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I.  INTRODUCTION 

A  need  exists  at  RADC  to  position  an  aircraft  with  respect  to  a  ground 
based  tracker.  Such  information  can  be  used  for  adjusting  and  subsequently 
calibrating  and  evaluating  the  tracking  performance  of  the  system.  A 
variety  of  operational  aircraft  will  be  used  which  .sharply  limit  the  equip¬ 
ment  that  may  be  conveniently  carried  by  the  aircraft.  Accordingly,  a 
passive  DME  approach  is  planned,  using  ground  based,  laser  equipped  units, 
reflecting  off  retroreflective  tape  affixed  to  the  skin  of  the  aircraft. 

The  objective  of  this  study  is  to  provide  an  analytical  tool  to 
assess  the  aircraft  positional  error  with  respect  to  the  tracker  position 
due  to  several  primary  error  sources.  This  report  discusses  the  concepts, 
theory  and  mathematical  formulations  which  lead  to  this  objective.  Sub¬ 
sequent  work  will  include  programming  of  these  results  for  graphical 
presentations  in  which  the  factors  of  error  magnitudes  and  network  geometry 
can  be  exercised. 

II.  TECHNICAL  DISCUSSION 

The  description  of  the  trilateration  network  is  made  with  reference  to 
Figure  II - 1 .  A  minimum  of  three  ground  based  DME  units  (i)  are  to  be  used 
to  passively  track  a  single  aircraft  (a/c) .  The  objective  of  this  study 
has  been  to  provide  the  conceptual  and  procedural  framework  for  assessing 
the  influences  of  the  various  error  sources  on  the  coordinates  of  location 
of  the  a/c  with  respect  to  the  tracker  at  any  time  (t) . 

The  definitions  of  symbols  and  terms  may  be  found  in  Appendix  B. 

A.  ERROR  SOURCES 

The  primary  errors  are  in  the  observations  of  range  by  the  DME 
units,  the  internal  survey  coordinate  errors,  timing  errors  and  the  bias 
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error  between  the  tracker  and  DME  local  survey  positions.  With  the 
exception  of  the  last  error,  which  is  assumed  to  be  composed  of  systematic 
biases (SX,  SY,  SZ) ,  the  errors  are  assumed  to  be  normally  distributed 
about  a  zero  mean  (random) .  The  character  of  the  errors  and  adopted 
symbolization  is  presented  in  Table  I I - 1 . 


ajc 


Figure  I I - 1  Trilateration  Network 
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ERROR  SOURCE 

RANDOM 

SYSTEMATIC 

Range  by  DME 

<5* 

Sr  =  o 

Survey  coord,  of  DME 

^  fifi  fill 

6XC  3  SY i  s  %Z'L  SO 

units 

Timing 

61 

St  s  o 

Survey  coord,  of 

=  6ie=0 

SXf  ,  SYe  ,  SZt 

tracker  w/r  DME 


TABLE  I I - 1  Error  Sources  and  Symbolization 


B.  MATHEMATICAL  MODELS 


The  models  required  for  the  error  analysis  relate  the  observed 
quantities  as  the  dependent  variables  to  the  unknown  and  given  quantities. 
Working  in  three-dimensional  rectangular  coordinates,  the  observations  of 
range  can  be  related  to  the  position  of  the  a/c  at  any  time  (t)  from  any 
DME  station  (i)  by: 


However,  all  the  DME  units  do  not  measure  range  simultaneously.  Conse¬ 
quently,  it  is  necessary  to  express  the  spatial  track  of  the  a/c  for  a 
short  period  in  terms  of  a  line  segment  in  space.  The  symmetric  form  of 
a  line  in  cartesian  space  is: 
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Xl  -  X2  Yi  -  Y2  Zi  -  Z2 
- a -  *  - 5 -  “  - c -  (2-2) 

where:  (1)  and  (2)  represent  two  points 

on  a  line  segment. 

We  now  assign  (X0,  Y0,  ZQ)  as  coordinates  of  the  a/c  at  about  the  midpoint 
of  the  line  segment  track.  At  any  other  time,  the  a/c  coordinates 
(Xt,  Yt,  Zt)  can  be  represented  as: 

Xt  =  X0  +  X  At  (2-3  a) 

Yt  =  Y0  +  Y  At  (2-3  b) 

Zt  =  Z0  +  Z  At  (2-3  c) 

where: 

•  •  • 

X,  Y,  Z  represent  the  time  rate  of 
change  of  a/c  along  the 
respective  axes 
At  =  time  at  ft)  -  tQ 
t0  *  the  time  when  the  a/c  is  at 
(X0,  Y0,  Z0) 

In  equation  (2-2),  the  points  can  be  replaced  by  the  origin  point  (XQ,  Y0, 
Z0)  and  any  other  point  at  time  (t) .  ie  (Xt,  Yt,  Zt) . 
i.e. : 

Xn  -  (Xp  +  XAt)  3  Y n  -  (Yo  ♦  YAt)  =  zn  -  (Zn  ♦  ZAt)  (2-4) 

a  b  "  c 

After  some  derivation  there  results: 

Xt  =  X0  +  %(Y  At) 

Yt  -  Yo  +  Jj(Z  At) 

zt  =  Zo  ♦  §(X  At) 
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(2-5a) 

(2-5b) 

(2-5c) 


Substituting  this  result  into  equation  (2-1): 
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f  (fit) =  Ttt  -  [  [Xl  “Xo 

[?.  -  -  s- ( 2>  n] 


('f-tU 


r 


[jc-Zc  -k(X*t)\T-0  M 


Equation  (2-6)  represents  the  mathematical  model  relating  observations  of 
range  by  DME  (i)  at  time  (t)  to  the  a/c  positional  parameters  of  the 
trilateration  network. 

To  facilitate  the  introduction  of  error  estimates  of  the  DME  station 
coordinates  in  a  local  survey  system,  it  is  convenient  to  treat  the  survey 
coordinates  as  if  they  were  directly  observed  quantities.  Their  weights 
can  then  be  used  in  the  adjustment  computation  to  introduce  their  contribu 
tion  to  the  system  error.  In  this  case  of  direct  observations  of  unknown 
parameters,  the  mathematical  models  are  simply: 

f(Xi)  =  -  Xi  =  0  (2-7a) 

obs  adj 

f (Yi)  =  Yi  -  Yi  =  0  (2-7b) 

obs  adj 

f(Zi)  =  Zi  -  Zi  =  0  (2-7C) 

obs  adj 

where : 


"obs"  pertains  to  the  observed  value 
"adj"  pertains  to  the  adjusted  value 


The  mathematical  models  which  relate  the  observations  of  DME  range 
and  local  survey  coordinates  to  the  parameters  defining  a/c  position  are 


the  Equations  (2-6)  and  (2-7). 

C.  OBSERVATION  EQUATIONS 

The  vehicle  through  which  this  analysis  of  errors  is  conducted  is  the 
adjustment  computation  using  least  squares.  For  this,  it  is  necessary  to 
arrange  the  mathematical  models  in  linear  form  with  respect  to  the  unknown 
parameters.  For  this,  the  model  is  represented  by  a  Taylor  Series  expansion 
neglecting  second  degree  and  higher  terms.  A  residual  is  appended  to  each 
observation  to  account  for  the  specific  observational  error.  These 
residuals  are  assumed  to  be  normally  distributed  about  a  zero  mean,  a 
condition  which  is  reasonably  well -attained  provided  the  mathematical 
models  represent  the  system.  In  this  case,  after  the  series  expansion  of 
the  models,  the  observation  equations  are:  [Merchant,  1974] 

V  +  BA  *  E  *  0  $0 

After  a  logical  grouping  of  parameters,  the  observation  equations 
pertaining  to  observations  of  DME  range  are: 

rs^PpPS-r 

V  +BA  *-BA  +BA  *  E  *  0  (a-?) 
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where:  (see  Appendix  B) 

r  _ 

V  =  vector  of  obs.  residuals  of  DME  range. 

b=  Wan  t) 

Hf(rlt)/5(Xo)Yo)Zo,^yo)Z.D 

A5  HXi  j  SY[  )5Z^]  for  all  ((.) 

a*  [ixo^Yoiszo)s>:o)sv'o)s±<Zl 

Ay  [sa^bjSc] 

r 

E  =  A  vector  of  f  )  evaluated  for  station  (i)  at  time 

(t)  using  the  observed  DME  range  and  current  estimates  of 
all  parameters  in  Equation  (2-6). 

The  observation  equations  pertaining  to  the  direct  observations  on  the 
local  survey  coordinates  of  DME  stations  are:  [Merchant,  1974] 

5  1  5 

V  -  A  *  E  =  0  (2-10) 

where : 

5 

V  =  vector  of  obs.  residuals  on  direct  survey  coordinate 

observations  (Xi,  Y^,  Zj.)  for  all  stations  (i) . 

S  _ 

E  =  a  vector  of  f  (X,  Y,  Z)  ^  evaluated  for  station  (i)  using 
original  obs.  and  current  estimates  of  all  parameters, 
[equations  (2-7)  ] 

Collecting  the  observation  equations  (2-9)  and  (2-10),  they  are 
expressed  as: 

V  +  B A  +  1=0  (2-11) 
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where : 


•£>  3 

O  O 
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r 

r 

V 

• 

5  V 

£ 

V  ‘ 

s 

LvJ 

[t. 

D.  NORMAL  EQUATIONS 

The  number  of  observation  equations  normally  exceeds  the  number  of 
unknown  parameters.  As  a  consequence,  Equations  (2-11)  have  no  unique 
solution.  By  imposing  on  these  equations  the  condition  that  the  sum  of  the 
weighted  squared  residuals  be  a  minimum  ("least  squares"),  a  solution  is 
uniquely  possible  provided  the  unknown  parameters  are  linearly  independent. 
This  condition  on  Equations  (2-11)  leads  to  the  "normal  equations": 

V  A  +  U  =  O  (a-/a) 

However,  since  we  are  interested  now  only  in  an  analysis  of  errors  as 
they  propagate  through  to  the  results,  we  need  not  solve  for  the  solution 
vector  (A).  We  need  only  form  and  invert  (N)  since: 


where : 


the  estimate  of  variance/ covariance  of  the  adjusted 


parameters . 

Accordingly,  we  are  interested  only  in  the  formation  and  inverse  of 
(N) ,  the  normal  coefficient  matrix.  It  can  be  shown  that:  [Merchant,  1974] 


N  *  F'w  B  (2-13) 

The  normal  coefficient  matrix  (N)  can  be  simplified  by  exploiting  the 
sparse  patterns  of  (B)  and  (W) .  The  result: 


(See  Appendix  A  for  derivations  of  elements  of  N) 

Provided  0,  a  unique  solution  for  the  system  of  18  equations  exists. 

The  computation  of  the  inverse  of  the  symmetrical  matrix  (N)  is  a  modest 
burden.  Accordingly,  no  further  attempts  are  made  to  further  exploit  the 
sparse  patterns  that  exist  in  (N)  for  purposes  of  computing  (N  **  ) . 

E.  RANDOM  POSITIONAL  ERRORS 

t 

By  extracting  the  array  (  2  )  within  (  2  )  related  to  a/c  position 
(ie  X0,  Y0,  Z0,  X0,  Y0,  Z0) ,  the  error  in  a/c  position  at  any  time  (t)  can 

p 

be  determined  by  propagation  of  ( 2  )  through  the  following  functions: 


Xf  -X0*Xo  At 

(a-/5a) 

X r  Yo^o 

(a  -/j  b) 

At 

(arl5c) 

F.  RANDOM  TIMING  ERRORS 

Errors  are  assumed  for  the  time  of  observation  (At)  as  well  as  for  the 
positional  parameters.  Therefore,  the  variance/ covariance  matrix  will  be: 

(A.  ~Ko) 

7X7 

G.  VARIANCE/COVARIANCE  OF  A/C  POSITION 

Finally,  the  estimated  variance/covariance  for  a/c  position  at  any  time 
(t)  due  to  errors  in  DME  observation,  time  observation  and  survey  (local) 
coordinates  can  be  computed  by  propagation  through  the  related  function. 

*G  i  G  (a -n) 
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where : 


* 


G  = 


2  * 


ie:  The  variance/covariance  matrix  for  the  spatial 

coordinates  of  the  a/c  at  the  sample  point  time  (t)  for 
error  sources  of: 

-  DME  range 

-  DME  spatial  coordinates 

-  time 

H.  RESULTING  RANDOM  AND  SYSTEMATIC  POSITIONAL  ERRORS 

The  foregoing  discussion  concerned  the  random  components  of  error 
producing  sources  in  the  trilateration  procedure.  The  primary  systematic 
error  sources  in  the  trilateration  procedure  are  the  fixed  or  bias  errors 
along  coordinate  axes  relating  the  local  survey  of  the  DME  gound  stations 
to  the  position  of  the  "tracker".  These  error  values  depend  on  the  means 
by  which  they  are  determined.  These  vary  with  each  application.  Since  the 
component  bias  errors  of  the  tracker  to  DME  distances  cannot  be  combined 
with  the  remaining  random  components  of  the  trilateration  procedure,  the 
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<rxt  vt  <Txtzt 


(Tv 


O' 


VfcZfc 

a 


3 *3 


dXt/a  (X0Y0ZoXoY0Zo  At) 
d\/d  (X0Y0Zo  XXZo  At) 
dZt/d(YXZoXXt,  At) 


total  system  error  for  a/c  position  at  any  time  can  be  expressed  as: 


Xt  error  * 

&Xt 

+ 

SXt 

(2-18a) 

Yt  error  ■ 

<ryt 

+ 

8  Y* 

(2-18b) 

Zt  error  = 

8Zt 

(2-18c) 

III.  CONCLUSION 

The  purpose  of  this  report  has  been  to  present  the  concept,  theory 
and  mathematical  basis  for  assessing  the  combined  errors  of  a  DME  ground 
based  tracking  system  on  the  relative  position  of  a  "target"  aircraft  with 
respect  to  a  ground  based  tracker. 

The  component  errors  of  the  system  that  have  been  included  are  those 
of  error  in  DME  passive  range  observations,  local  geodetic  survey  coordinates 
of  the  ground  based  DME  units  and  system  timing  errors.  The  result  is  in 
terms  of  an  algebraic  expression  of  the  variance/covariance  of  the  aircraft 
with  respect  to  the  tracker's  position.  In  addition,  provision  was  made 
for  the  introduction  of  the  geodetic  bias  error  which  exists  between  the 
tracker  and  the  local  DME  ground  survey. 

Subsequent  work  to  meet  the  overall  objective  should  include,  in  order: 

-  analysis  of  the  design  matrix  (B)  to  assure  linear 
independence 

-  programming  of  the  Equations  (2-14,  and  2-17) 

-  generating  plots  of  positional  errors  exercising  the  geometry 
of  the  network  and  the  magnitude  of  contributing  error  values 
which  are  representative  of  anticipated  applications. 
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APPENDIX  A 


ELEMENTS  OF  THE  DESIGN  AND  WEIGHT  MATRICES 


-  for  local  DME  survey  coordinates 

-  further  breakdown: 


where  #  DME  sta  i=  1  to  n 
for  ^  3 

ii  *  [if  (ra)  /d(xtYtzj] 

m 


where  j  =  1  to 

for  DME  obs.  from 
sta.  (i). 
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for  parameters  of  a/c  position 
further  breakdown: 


ir&f  (rL])/^UxZoX%i)} 

j 

le:  Si 

3i  (i>)  '('zUr^/zz,) 
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for  parameters  describing  spatial  line  segment 
further  breakdown: 


r 

W  -  Weight  matrix  for  DME  range  observations: 


Ideally  p  f 

U  -  2" 


However,  since  correlations  between  observations  of  range  should  be  low, 
provided  careful  modeling  has  been  accomplished,  the  weight  matrix  can  be 
taken  as  a  diagonal  matrix  whose  elements  are  the  inverse  of  the  standard 


error  (standard  deviation),  squared  in  the  units  in  which  the  observation 


where : 
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I  -  weight  matrix  for  local  survey  coordinates  of  the  DME  ground  station 
locations 

$  A*l 

Ideally,  W  *  2*  ;  i.e.,  the  inverse  of  the  full  variance/covariance 

matrix  which  relate  all  station  coordinates.  However,  a  valid  estimate 
of  (£)  is  rarely  available.  Consequently,  the  correlation  between 
coordinate  values  of  survey  points  is  normally  neglected.  Accordingly; 


APPENDIX  B 


DEFINITIONS  OF  TERMS 

aircraft 

standard  deviation  of  normally  distributed  quantity 

about  a  zero  mean 

variance 

covariance  of  x  with  y 
systematic  bias  error 

coordinates  in  a  three  dimensional,  rectangular, 
right  handed  system 
range  time 

constants  appearing  in  the  symmetric  form  of  an 

equation  expressing  line  in  space 

position  of  a/c  at  time  (t) 

position  of  a/c  at  time  (to) 

variance/covariance  matrix 

coordinates  of  DME  station  (i) 
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selected  acquisition  programs  in  support  of  Command,  Control 
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